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**PREFACE** 


This  report  presents  the  results  of  an  independent  investigation 
into  the  accuracy  of  horizontal  sextant  angles  as  a  method  of  establish¬ 
ing  geographic  position  conducted  by  the  authors  from  October  1973  to 
March  1974.  It  includes  a  description  of  the  system,  the  development  of 
a  method  for  determining  positioning  accuracy  and  an  example  problem  with 
a  step-by-step  procedure  for  utilizing  the  developed  method.  Supporting 
theory,  analysis  and  proofs  are  included  wherever  possible  for  the  con¬ 
venience  of  the  reader. 
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SUMMARY 


The  text  is  preceded  by  a  brief  introduction  which  outlines  the 
history  of  horizontal  sextant  angles  as  a  positioning  method,  the  reason 
for  the  recent  interest  in  defining  the  accuracy  of  the  system  and  the 
general  direction  of  the  analysis  which  follows.  The  description  of  the 
method  of  horizontal  sextant  angles  as  a  means  of  establishing  geographic 
position,  including  the  theoretical  basis  for  its  use  and  an  illustra¬ 
tion  of  a  particular  method  of  utilization,  is  contained  in  Chapter  1. 

A  general  description  of  the  instruments  used  is  included  also.  Chapter  1 
ends  with  the  development  of  the  error  model  describing  each  of  the 
component  dependent  and  independent  angular  errors  which  affect  an  exact 
angular  measurement.  Errors  in  the  use  of  both  the  three-arm  protractor 
and  the  sextant  are  considered.  The  error  model  is  used  in  Chapter  2  to 
develop  the  basic  method  for  determining  positioning  accuracy. 

The  development  of  the  basic  method  for  determining  positioning 
accuracy  is  begun  with  an  explanation  of  the  necessary  mathematical 
considerations  and  a  brief  look  at  the  applicable  probability  theory. 

The  steps  in  the  actual  development  include: 

1.  The  definition  of  the  individual  LOP  (line  of  position) 
variances  and  the  combined  LOP  covariance  from  the  component 
angular  errors. 

2.  The  conversion  of  the  variances  and  the  covariance  to 
distance  units. 

3.  The  elimination  of  the  dependence  of  the  LOP  errors  using 
Fraser's  method  (ref.  3). 


Preceding  page  blank 
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4.  The  determination  of  the  radius  of  the  circular  probability 
density  contour  within  which  the  actual  position  can  be 
expected  to  lie. 

The  Geometric  Multiplication  of  Precision  (GMOP)  is  derived  and 
analyzed  in  Chapter  3.  The  GMOP  is  defined  as  a  dimensionless  quantity 
ft  per  unit  baseline  length,  where  K  is  the  conversion  factor  used  to 
obtain  the  LOP  variances  and  covariance  in  distance  units.  (The  con¬ 
version  is  from  radians  to  distance  units.)  The  similarity  between 
the  characteristics  of  the  gradient  of  ji  is  defined  as  the  exact 

angular  measurement  of  a  horizontal  sextant  angle)  and  the  desired 
traits  of  the  GMOP  sets  the  direction  of  the  development  toward  finding 
(MOP  as  a  function  of  grad ^3  .  The  equation  of^  in  the  polar  coordi¬ 
nate  system  is  found,  and  the  GMOP  is  derived  as  the  reciprocal  of  the 
magnitude  of  grad j$.  The  GMOP  is  analyzed  thoroughly  to  investigate  its 
physical  correlation  and  its  functional  behavior. 

Chapter  4  contains  the  derivation  of  a  simple  method  for  determin¬ 
ing  tne  IX)P  crossing  angle  oj.  .  This  quantity  ir  required  for  the 
implementation  of  the  basic  method  described  in  Chapter  2.  It  is  shown 
that  is  easily  obtainable  from  the  measurable  quantities  which  are 
required  for  the  calculation  of  the  GMOP. 

Chapter  5  presents  an  example  problem  for  the  convenience  of  the 
reader  with  a  step-by-step  procedure  for  applying  the  total  method  to 
a  practical  situation.  The  lack  of  suitable  published  human  factors 
data  for  the  three-arm  protractor  and  the  sextant  forces  certain 
assumptions  concerning  the  statistical  distributions  of  the  component 
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angular  errors.  All  assumptions  are  explained  before  the  problem 
solution  Is  presented. 

Appendix  A  offers  the  methods  of  obtaining  the  component  angular 
error  distributions  for  the  human  factors  errors  In  using  the  three-arm 
protractor  fiom  measurable  distance  errors  in  placing  the  pivot  point 
and  the  arms.  These  methods  can  be  applied  to  suitable  experimental 
data  to  obtain  the  component  angular  error  distributions  from  which 
the  probabilistic  variables  necessary  for  the  implementation  of  the 
method  for  determining  positioning  accuracy  are  derived. 

Appendix  B  contains  the  proof  of  the  acceptability  of  the  LOP-error 
model  of  fig.  2-1.  Specifically,  the  contention  that,  in  the  neighbor¬ 
hood  of  the  observers'  position,  the  LOP's  approximate  straight  lines  and 
the  grad curves  along  which  the  error  is  measured  are  essentially 
straight  lines  perpendicular  to  the  exact  and  the  error  LOP's  is  proven 
to  be  reasonable.  Appendices  C  and  D  show  some  intermediate  steps  for 
derivations  contained  in  Chapter  3. 
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INTRODUCTION 


The  advantages  in  convenience  and  precision  of  using  horizontal 
angles  for  finding  a  position  on  land  was  recognized  certainly  before 
the  seventeenth  century.  The  application  of  the  method  to  hydrographic 
surveys  and  shipboard  navigation  was  delayed  considerably  due  to  the 
lack  of  a  suitable  angle  measuring  Instrument.  The  invention  of  the 
octant  by  John  Hadley  in  1731  provided  such  an  instrument  for  marine 
use,  and  the  Rev.  John  Mlchell  suggested  the  use  of  the  octant  for 

Q 

triangulation  afloat  in  1765. 

Since  the  time  of  Hadley  and  Mlchell,  the  application  of  horizontal 
sextant  angles  to  marine  positioning  problems  which  require  great 
precision  has  increased  markedly.  One  specific  application  for  which 
the  method  is  suited  well  is  buoy  anchor  positioning.  The  U.  S.  Coast 
Guard  maintains  more  than  20,000  floating  aids  to  navigation  which  mark 
navigable  channels  and  hazards  to  navigation  in  the  waters  of  the 
United  States.  The  Coast  Guard  Aids  to  Navigation  Manual  (CG-222) 
prescribes  horizontal  sextant  angles  as  the  primary  buoy  anchor  position¬ 
ing  method  to  be  used  by  U.  S.  Coast  Guard  buoy  tenders.  Interest  in 
the  accuracy  of  this  positioning  method  has  been  spurred  within  the 
Coast  Guard  by  the  growing  population  of  very  large  merchant  vessels 
whose  requirements  for  accurate  channel  marking  are  emphasized  by  the 
extremely  severe  consequences  of  grounding,  e.g.,  cargo  loss,  ecological 
catastrophe  and  vessel  damage. 

A  method  for  determining  the  accuracy  which  can  be  expected  for  a 
particular  horizontal  sextant  angle  position  fix  is  derived  and  presented 
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in  this  report.  Chapter  1  offers  the  theoretical  basis  for  using 
horizontal  sextant  angles  as  a  positioning  method  and  the  development 
of  the  system  model  for  a  three-object,  two-llne-of-posltion  fix. 

The  basic  method  Is  developed  from  a  random  error  analysis  in  Chapter  2. 
Chapters  3  and  4  and  Appendix  A  contain  supporting  derivations  which 
provide  necessary  inputs  to  the  basic  method  of  Chapter  2.  An  example 
problem  is  presented  in  Chapter  5  for  the  convenience  of  the  reader  to 
demonstrate  the  application  of  the  total  method  to  a  practical  situation, 
including  a  step-by-step  procedure  which  facilitates  the  calculations. 


i 


-18- 


1.  HORIZONTAL  SEXTANT  ANGLES  AS  A  POSITIONING  METHOD 
1.1  Theoretical  Basis 

The  utility  of  horizontal  sextant  a igles  as  a  means  of  establishing 
geographic  position  (called  position  fixing)  is  based  upon  the  following 
geometric  principles: 

a.  Any  three  points  in  space  are  coplanar  and  define  a  circle 
whose  orientation  and  dimensions  depend  upon  the  relative 
positions  of  the  three  points. 

b.  Consider  a  circle  with  two  fixed  points  on  its  circumference 
which  are  connected  by  a  chord.  The  angle  formed  by  the  line 
segments  constructed  from  each  of  the  two  fixed  points  to  any 
point  P  on  the  circumference  is  constant  for  all  points  P 

in  the  same  segment  of  the  circle.  (A  segment  of  a  circle  is 

that  portion  of  the  circle  which  is  bounded  by  ar.  arc  and 
1 

its  chord.  ) 

The  fact  that  any  three  points  in  space  define  a  plane  is  basic  and 
will  not  be  justified  further  herein.  The  proof  of  the  second  part  of 
principle  "a"  is  founded  upon  the  equidistance  of  all  the  points  on  the 
circumference  of  a  circle  from  its  center.  Fig  1-1  shows  a  typical 
triangle  OPQ.  The  perpendicular  bisectors  of  OQ  and  OP,  labeled  WX  and 
YZ  respectively,  are  shown  also.  The  end  points  of  the  bisected  segments 
are  equidistant  from  any  particular  point  on  the  perpendicular  bisector. 
Therefore,  0  and  Q  are  equidistant  from  any  particular  point  on  WX;  and 
0  and  P  are  equidistant  from  any  particular  point  on  YZ.  It  follows  that 


FIG.  1-1  TYPICAL  TRIANGLE  OPQ  SHOWING 
CIRCUMSCRIBING  CIRCLE 


P 


FIG.  1-2  DIFFERENT  TRIANGLES  WITH  COMMON 

SIDE  OQ  AND  THE  SAME  CIRCUMSCRIBING 
CIRCLE 
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0,  P  and  Q  are  equidistant  from  the  point  of  intersection  of  WX  and  YZ, 
labeled  C.  Consequently,  0,  P  and  Q  are  on  the  circumference  of  a 
circle  whose  center  is  C.  Further,  since  WX  and  YZ  are  straight  lines, 
they  intersect  at  only  one  point;  so  only  one  circle  is  defined  by 
0,  P  and  Q. 

The  contention  of  principle  "b"  is  illustrated  in  fig  1-2.  If 
0  and  Q  are  fixed  on  the  circumference  of  the  circle,  the  value  of  the 
angle  formed  by  line  segments  OP  and  QP  is  constant  for  all  points  P 
in  the  same  segment  of  the  circle,  i.e.,  ".  Fig  l-3a 

depicts  the  situation  when  the  lengths  of  OP  and  QP  are  equal  (OPQ  is 
an  isosceles  triangle).  The  isosceles  triangle  OCQ,  whose  sides  are 
the  baseline  0Q  and  the  radii  CO  and  CQ,  is  shown  also.  (The  standard 
symbols  A  for  triangle  and  Z  for  angle  will  be  used  henceforth.) 

Now,  the  line  segment  PS  is  constructed  as  the  perpendicular  bisector  of 
0Q,  which  is  the  base  of  AOPQ  and  a  chord  of  the  circle.  Consequently, 

PS  passes  through  C  and  bisects  A  OPQ.  Therefore,  ZOPC  =  ZOPS  =  ^9/2; 
and,  since  AOPC  is  also  isosceles,  Z  POC  =  ZOPC.  It  follows  that 
Z  OCP  =  TT  -  £  and  /.DCS  -f3  .  This  implies  that  Z  COQ  =  3L. 

Fig  l-3b  shows  another  possible  configuration  of  the  inscribed  triangle. 
Again,  0  and  Q  are  the  fixed  points  on  the  circumference  of  the  same 
circle;  but  this  time  the  triangle  is  defined  by  0,  P'  and  Q.  The 

triangles  OCQ  of  figs  l-3a  and  l-3b  are  identical.  Since  Z  COQ  =  ^  -A 

2  r 

from  fig  l-3b. 


Z  P'OC  ■  0  -  Z  COQ  -  6  +  -  -3l_  Z0P‘C 


(1.1) 
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Also,  ZOQP'  ■  'TT  -  0  -  ^  •  But  ZCQO-  Z  COQ,  so  that 

ZCQP*  -  ZOQP'  -  ZCQO  -  3L  -  0  -  0  '  +  yg  -  ^CP'Q  (1.2) 

By  definition, 

-  £  OP'C  +  ^  CP'Q  (1.3) 

Substituting  eqs.  1.1  and  1.2  into  eq.  1.3, 


V 


TT 

2 


JL  - 
2 


or  =  yg 


This  result  can  be  established  in  the  same  manner  for  all  other  points  P 
in  the  same  segment  of  the  circle. 

The  important  implication  of  the  preceding  developments  is  that,  for 
any  set  of  fixed  points  0  and  Q,  a  particular  value  of  p  defines  a 
unique  segment  of  a  circle  whose  arc  is  a  contour  of  constant  p  .  There¬ 
fore,  if  the  points  0  and  Q  are  fixed,  charted  objects,  the  angle  ^  observed 
from  any  position  P  defines  a  circular  line  of  position  (LOP)  which  passes 
through  P.  Further,  if  two  or  more  such  angles  are  measured  from  the  same 
position  P  using  different  objects  0  and  Q,  the  intersection  point  of  the 

I 

arcs  of  the  circle  segments  circumscribing  AOPQ,  AO'PQ',  AO"PQ"> 
etc. ,  will  fix  the  position  of  P  within  the  limits  of  accuracy  of  the 
system.  This  is  the  theoretical  basis  for  position  fixing  using  the  method 
of  horizontal  sextant  angles.  ' 


1.2  Practical  Application 

The  method  of  horizontal  sextant  angles  is  used  as  the  primary 
buoy  anchor  positioning  method  aboard  most  U.  S.  Coast  Guard  buoy  tenders. 


In  this  application,  the  position  fix  is  obtained  from  two  LOP's  de¬ 
rived  from  three  fixed  objects.  (The  reliability  of  the  fix  often  is 
assessed  by  comparison  with  a  yyro  bearing  fix.)  The  case  for  the  three 
object,  two-LOP  fix  versus  fixes  from  other  possible  combinations  of 
objects  and  LOP's  will  not  be  argued  in  this  report  in  order  to  preserve 
the  continuity  of  the  discussion.  To  align  the  discussion  with  the 
common  practice  for  positioning  buoy  anchors,  only  the  three-object, 
two-LOP  fix  will  be  considered  hereafter.  The  two  general  geometric 
situations  which  may  result  when  two  LOP's  are  obtained  from  three  ob¬ 
jects  are  shown  in  figs  1-4. 

The  example  of  utilizing  the  method  of  horizontal  sextant  angles 
which  will  be  described  is  called  the  preset  angle  method.  (The  choice 
of  this  particular  method  for  the  example  does  not  degrade  the  general¬ 
ity  of  subsequent  developments.)  The  mechanics  of  applying  this  method 
to  the  buoy  anchor  positioning  problem  begins  sometime  before  the  actual 
operation.  Files  and  notes  of  previous  positioning  operations  for  the 
same  buoy  anchor  are  researched  to  determine  the  three  fixed  objects 
(0,  Q  and  Q')  which  are  best  situated  for  the  particular  intended  anchor 
position  P.  Generally,  alternate  sets  of  objects  are  selected,  if 
available,  to  permit  some  flexibility  if  the  primary  objects  are  found 
to  be  unsuitable  at  the  time  of  the  positioning  operation.  Considera¬ 
tions  in  selecting  fixed  objects  include: 

a.  visibility  from  the  intended  buoy  anchor  position, 

b.  definition  of  vertical  centerline,  i.e.,  a  spire  is  more 
desirable  than  a  tangent  to  a  point  of  land, 
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c.  suitability  of  the  angle  of  intersection  (crossing  angle) 
of  the  LOP's,  and 

d.  the  observed  angle  must  not  exceed  the  maximum  value  which 
can  be  measured  by  the  sextant  (126°  for  the  newest 
sextants) . 

Relatively  elementary  instruments  are  used  to  measure  the  sextant 
angles  to  be  observed  and  the  LOP  crossing  angle.  The  three-arm  pro¬ 
tractor  is  used  to  measure  the  angles  from  a  standard  National  Ocean 
Survey  (NOS)  chart.  As  its  name  implies,  this  is  a  three-armed  instru¬ 
ment.  Its  arms  may  be  rotated  about  a  common  pivrt  point,  and  the 
angles  between  each  outer  arm  and  the  center  arm  are  read  directly  from 
a  fixed  scale.  The  vernier  scale  offers  marked  precision  to  one  or 
two  minutes  depending  upon  the  protractor  used.  There  is  a  small  hole 
at  the  pivot  point  which  is  centered  over  the  charted  buoy  anchor  posi¬ 
tion.  Then,  the  arms  are  aligned  through  the  charted  position  of  each 
of  the  three  fixed  objects.  The  measured  angles  are  recorded.  In 
general,  the  largest  scale  chart  available  upon  which  the  necessary  ob¬ 
jects  are  shown  is  used;  and  the  angular  measurements  are  checked  by 
several  competent  personnel.  The  LOP  crossing  angle  is  observed  usually 
by  constructing  the  circular  LOP's  on  the  chart  with  a  compass. 

As  the  vessel  proceeds  to  the  area  in  which  the  buoy  anchor  will 
be  positioned,  each  of  che  horizontal  sextant  angles  to  be  observed  is 
set  on  a  sextant.  Two  competent  personnel,  each  with  one  of  the  preset 
sextants,  assume  a  prescribed  position  on  the  vessel  which  is  selected 


for  its  proximity  to  the  point  from  which  the  anchor  will  be  released 
and  for  its  commanding  view  of  the  fixed  objects.  The  vessel  is 
maneuvered  until  the  vertical  centerlines  of  the  objects  used  to  obtain 
each  LOP  appear  superimposed  when  sighting  through  the  appropriate  sex¬ 
tant.  When  this  occurs,  the  anchor  is  released;  and  the  positioning 
operation  is  complete. 

1  -3  Sources  of  System  Error 

The  sources  of  error  in  the  system  as  described  in  this  chapter 
can  be  classified  under  three  general  headings.  These  are  plotting 
error,  primary  instrument  error  and  supplementary  error.  The  specific 
components  of  each  general  heading  are  the  errors  which  must  be  con¬ 
sidered  for  the  analysis  of  the  positioning  accuracy  of  the  system. 
These  include  the  following: 

(1)  Plotting  error 

(a)  charting  error  of  the  positions  of  objects 

(b)  instrument  error  in  angular  measurement  with  the 
three-arm  protractor 

(c)  human  error  in  angular  measurement  with  the  three- 
arm  protractor 

(2)  Primary  instrument  error 

(d)  instrument  error  in  angular  measurement  with 
the  sextant 

(e)  human  error  in  angular  measurement  with  the 
sextant 

(3)  Supplementary  error 

(f)  error  introduced  by  the  displacement  of  the 


sextant  angle  observers  from  the  position  of 
buoy  anchor 

(g)  error  caused  by  the  effect  of  current  on  the 

anchor  after  it  is  released  and/or  the  effect  of 
bottom  topography  on  the  position  of  the  anchor 
after  it  contacts  the  bottom. 


The  supplementary  error  will  be  neglected  in  the  error  analysis,  since 
its  components  are  not  truly  random. 

First  consider  the  combined  effect  of  (a)  and  (c).  This  combined 
error  is  the  source  of  dependent  error  between  the  two  angles  due  to  the 
use  of  the  three-arm  protractor.  Fig  1-5  shows  a  typical  set  of  fixed 
objects  0,  Q  and  Q1  which  are  to  be  used  to  determine  the  position  P. 

Line  segments  PO,  PQ  and  PQ'  represent  the  arms  of  the  three-arm  protrac¬ 
tor  used  to  measure  the  sextant  angles  to  be  observed.  ^  and  ^  are 
the  exact  angular  measurements  of  £0PQ  and  £  QPQ'  respectively.  The 
error  in  positioning  the  pivot  point  of  the  protractor  directly  over  P 


Pi 


as  it  affects  the  angular  measurements  p  ^  and  ^ ^  is  designated  £ 

and  £  .  Failure  to  position  the  arms  through  the  exact  positions 

p2 

0,  Q  and  Q1  creates  the  errors  £  ,  £  and  £  ,  in  the  magnitude  of 

o  q  q 

the  angle  measured  from  a  reference  direction  (REF)  to  arms  PO,  PQ  and 
PQ' .  The  angle  from  REF  to  each  arm,  indicated  in  fig  1-5,  can  be  sub¬ 
tracted  to  give  ^ OPQ  and  ^ QPQ' ,  which  define  LOP-1  and  LOP-2  respec¬ 
tively,  as 


£0PQ-  C  (REF-to-PO)  +  £  -  ^(PSF-to-PQ)  -  £  +  £ 

°  q  Pi 


(1.4) 


and 


£QPQ'  .  ^  (REF-to-PQ)  +  £  -  ^  (REF-to-PQ')  -  £  ,  +  £ 


p2 


-29- 


The  signs  of  the  errors  may  be  either  positive  or  negative.  Positive 

errors  £^,  and  tend  to  increase  the  magnitude  of  the  angle 

between  REF  and  the  particular  arm  while  positive  errors  £  and  £, 

P1  p2 

increase  the  total  angle.  The  dependence  of  the  error  terms,  shown  for 
both  angles  in  eqs.  1.4,  upon  each  other  is  apparent  from  the  appearance 
of  the  £q  and  £p  terms  in  both  equations. 

The  primary  instrument  error  and  (b),  the  three-arm  protractor 
instrument  error,  are  individually  independent  of  all  other  errors. 

Their  combined  effect  on  and  ?2  is  designated  £r  and 
respectively.  Adding  these  errors  and  noting  that 


^  (REF-to-PO)  -  £  (REF-to-PQ) 

and  ^  (REF-to-PQ)  -  ^(REF-to-PQ'), 

eqs  1.4  become 

<£0rq  -  fi J  +  C0  -  £  +  €  +  € 

r  ’  P1  (1.5) 

and  ^QPQ'  -  /S2  +  £,  -  +  £p2  +  £s 

The  eqs.  1.5  are  the  expressions  for  the  effective  horizontal  sextant 
angles  used  for  the  positioning  operation  described  in  Section  1.2.  The 
remainder  of  this  presentation  is  dedicated  to  the  determination  of 
positioning  accuracy  from  the  expected  angular  errors  £  . 
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2.  POSITIONING  ACCURACY 

2.1  Preliminary  Considerations 

The  ultimate  goal  of  this  report  is  to  develop  a  method  for  determin¬ 
ing  the  accuracy  of  a  position  fix  obtained  by  ho~izontal  sextant  angles. 
The  basic  method  is  derived  and  presented  in  section  2.3.  The  assumptions 
which  are  made  to  permit  a  reasonable  mathematical  development  of  accuracy 
probability  are  the  following: 

a.  the  combined  system  random  error  is  normally  distributed; 

b.  any  bias  error  has  been  removed  from  the  system  so  that 
the  mean  of  the  random  error  is  zero; 

and,in  the  neighborhood  of  the  observers'  position  P, 

c.  the  LOP's  are  coplanar; 

d.  the  exact  and  error  LOP's  approximate  straight  lines  and 
the  angular  errors  €  can  be  transformed  into  distance 
errors  measured  along  a  straight  line  which  is  perpendicu¬ 
lar  to  the  exact  LOP  and  the  error  LOP;  and 

e.  the  error  LOP's  are  parallel  to  the  exact  LOP's. 

The  statistical  descriptions  of  the  individual  angular  errors 

are  derived  from  analyses  of  suitable  samples  of  random  error  data.  (The 

determination  of  the  angular  error  distributions  of  those  errors  which 

result  from  the  use  of  the  three-arm  protractor,  £  Q,  £q,  £qi>  £.p 

and  £.  ,  from  measurable  distance  errors  in  placing  the  arms  and  the 

p2 

pivot  point  is  discussed  in  Appendix  A.)  Generally  speaking,  random 
physical  phenomena  approximate  normal  distributions.  Further,  the  sum  of 
random  variables  with  the  same  distributions  which  are  other  than  normal 
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usual  ly  will  approach  a  normal  distribution.  Consequently,  the  assump¬ 
tion  "a"  of  a  normally  distributed  random  error  for  the  horizontal  sex¬ 
tant  angle  problem  is  considered  acceptable.  The  distribution  of  some 
random  error  variables. is  offset  from  a  zero  mean  value  by  a  constant 
value  called  the  bias  error.  In  order  to  concentrate  the  discussion  on 
the  analysis  of  the  system  random  error,  it  is  assumed  by  "b"  that  the 
bias  error  has  been  removed  from  the  system. 

The  geometric  model  used  for  this  analysis  is  shown  in  fig.  2-1. 

This  model  is  identical  to  that  used  by  Burt  et  al  in  ref.  2.  The  fig¬ 
ure  depicts  an  enlarged  view  of  a  small  area  in  the  neighborhood  of  the 
observers'  position  P.  The  effect  of  a  one  standard  deviation  (one  -  <7"  ) 
error  in  each  LOP  is  shown.  The  combination  of  the 

one  -  <J"  error  LOP's  for  the  assumed  normally  distributed  random  error 
defines  the  elliptical  probability  density  contour  of  a  bivariate  normal 
distribution.  The  ellipse  shown  in  fig.  2-1  might  represent  the  757. 
probability  density  contour,  i.e.,  there  is  a  75%  probability  (three  chances 
in  four)  that  any  particular  position  fix  will  be  within  its  boundaries. 
Assumptions  "c" ,  "d"  and  "e"  are  necessary  for  the  use  of  this  model. 

The  coplanarity  of  the  LOP's,  assumption  "c",  is  considered  accept¬ 
able  if  the  fixed  objects  are  clearly  visible  from  the  observers  position 
and  if  the  fixed  objects  and  the  observers  position  are  at  the  same 
altitude.  The  first  requirement  is  implicit  for  the  application  of  the 
method  of  horizontal  sextant  angles.  The  latter  is  usually  reasonable 
for  buoy  anchor  positioning  problems.  The  linearity  of  the  LOP's  and 
the  distance  error  lines  assumed  by  "d"  is  made  temporarily  to  preserve 


the  continuity  of  the  discussion.  Proof  of  the  acceptability  of  this 
assumption  is  presented  in  Appendix  B.  Some  of  the  developments  of 
Chapter  3  are  used  in  the  proof.  Fig.  2-2  shows  that  the  assumption  "e" 
of  parallel  exact  and  error  LOP's  is  relatively  poor  near  the  points  0 
and  Q.  For  the  relatively  small  expected  LOP  error,  the  LOP's  will  be 
essentially  parallel  everywhere  except  in  small  areas  in  the  neighbor¬ 
hoods  of  0  and  Q.  Consequently,  this  assumption,  with  its  recognized 
shortcomings  near  the  fixed  objects,  is  considered  acceptable  for  the 
purposes  of  this  anel/sis. 

2.2  The  Probability  Theory 

The  probability  density  function  for  a  normally  distributed  random 
error  variable  £  with  a  zero  mean  value  is  given  by 


p(  e  )  - 


^2=rF<r 


e 


1st 

2  <r*- 


(2.1) 


where  (J~  is  the  standard  deviation  of  the  distribution  and  is 


the  variance  of  the  distribution. 

The  probability  that  a  particular  value  of  £  is  within  the  limits  £ 
and  is  determined  by  integrating  eq.  2.1  from  £a  to  In 

equation  form, 

p  (  ca  <  e  ^  eb>  =  J£a  p(  £) (2-2) 

Tabular  data  is  available  for  the  evaluation  of  eq.  2.2  for  any  particular 

values  of  £  ,  £,  and  (T ' . 

a  d 
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The  general  form  of  the  joint  probability  density  function  of  two 
normally  distributed  random  error  variables  £ ^  and  £  ,  both  with 
zero  mean  values,  is  the  bivariate  normal 


1  e  2(1  -  f*> 

2-rrcjjo^i  -  f 


where  ^  is  the  correlation  coefficient  (a  measure  of  the  dependence) 

of  CL  upon  £  • 

The  probability  that  a  particular  jalue  of  the  cabined  error  is  within 
certain  prescribed  limits  of  ^  and  can  ^e  determined  by  the  double 

integration  of  eq.  2.3  between  these  limits,  i.e., 
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The  dependent  nature  of  each  total  LOP  error  upon  the  other  tor  the 

three-object,  two-LOP  method  of  horizontal  sextant  angles  was  discussed 

in  section  1.3.  This  dependence  requires  the  more  complicated  orobability 

density  function  of  eq.-2.3.  Fraser  (ret.  3)  gives  a  method  by  which 

the  dependent  errors  £  ^  and  £^  can  be  transformed  into  orthogonal, 

independent  errors,  x  and  y,  whose  probability  density  function  is  in 

the  form  of  eq.  2.5.  Fig.  2-3  (fig.  12  of  ref.  2)  can  be  entered,  then, 

directly  with  a  prescribed  function  of  QT  and  to  obtain  the  radius  of 

x  y 

an  equivalent  circular  probability  contour  for  various  accuracy  probabili¬ 
ties.  (The  elliptical  probability  density  contours  of  the  bivariate 
normal  are  converted  to  circular  contours  in  order  to  quantify  more  de¬ 
scriptively  the  accuracy  to  be  expected  for  a  particular  probability. 

The  actual  dimensions  of  the  error  ellipse  are  of  interest  only  when  a 
specific  situation  requires  the  determination  of  along-track  and/or  cross¬ 
track  error.) 

A  suitable  modification  of  Fraser's  method  for  application  to  hori¬ 
zontal  sextant  angles  and  the  subsequent  use  of  fig.  2-3  will  be  combined 
to  form  a  method  by  which  the  positioning  accuracy  of  the  system  can  be 
determined.  The  first  task  is  to  define  the  variances  (<r2)  of  the  indi¬ 
vidual  LOP  error  and  the  correlation  coefficient  (^)  of  the  combined 
LOP  error. 

2.3  The  Basic  Method 

The  analysis  begins  with  the  definition  of  the  equations  for  each 


LOP  from  eqs.  1.5  as 
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lop  - 1  -  z  opq  -  ^  +  e0  -  £q  +  cp  +  £r 


(2.6) 


and  LOP  -  2  -  Z  QPQ'  »  fS2  +  Gq 


+  € 


s 


As  discussed  in  Section  1.3,  the  error  components  combined  as  £r  and  £g 
are  individually  independent  of  all  other  errors.  Further,  the  errors  in 
positioning  the  outer  arms  of  the  three-arm  protractor,  £q  and  £qt , 
independently  affect  LOP  -  1  and  LOP  -  2  respectively.  These  independent 
errors  are  combined  to  facilitate  the  mathematical  manipulations  as 

£t  -  CQ  +  £r 

and  £h«-£qi  +  £s 

Now,  the  expressions  for  the  total  error  affecting  each  LOP  can  be  written 
as 

^1  “  ^LOP  -  1  ^i  "  ^-q  +  ^■pi  (2.7a) 

and  C2  =  Clop  .  2  =  £h+  £q+  £p2  (2.7b) 


By  definition,  the  variance  of  £  is  equal  to  the  expected  value  of  the 
square  of  £  .  In  equation  form 

(T2*e[€2]  (2.8) 

Applying  eq.  2.8  to  eq.  2.7a, 

<ri2-*[<ei-  V  Vi 
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or 

or 


l"  *  E  ^-q  +  +  £, 


-  £  £ 

q  q  p. 


C4  C  -  £ 


Pi  q 


S +  €fi  ] 


(2.9) 


Since  E^X  +  Y  +  Z  + _ ]  “  E[  X]  +  E[Y]  +  e[  Z  +  ...]  ,  eq.  2.8  can 

be  applied  to  eq.  2.9  to  yield 

2  (2.10) 


«7 


^  ^  +  S  '2E[  £‘V  €«  S'  £i  SJ 


The  last  three  terms  in  eq.  2.10  are  the  covariance  (j/jC)  terms  which  express 
the  dependence  of  the  random  variables .  However,  as  stated  previously,  the 
error  is  individually  independent  of  all  other  errors,  so  that 
-  E£  £,t  Cq  J  -  0  and^6tip  =  E  £  C.^  £p  j  =  °.  Also,  assuming 
that  the  error  in  aligning  arm  PQ  directly  through  Q  is  independent  of  that 
in  positioning  the  pivot  point  directly  over  point  P,  W  B  e|*  £  £  ■  0 

^  qpi  [  q  piJ 

and  all  the  covariance  terms  vanish.  Eq.  2.10  becomes 


or2+  crz+  cr 


(2.11) 


A  similar  sequence  of  operations  on  eq.  2.7b  yields 


or2-  cr2  +  cr2+  cr 


(2.12) 


The 


correlation  coefficient  ^  of  £i  and  ^  *-s  defined  as 


r  a^cr[ 


(2.13) 
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where  *  E^  €2J  is  the  covariance  of  £j  and  ^  .  (2.14) 

Substituting  eqs.  2.7  into  the  expression  for  the  covariance,  eq.  2.14, 


>*12*  E[  <  €i  '  s +  S)<£h+  £’+  V] 

or 

/*12  [  1  h  *  ^  1  P2  *  h  *  1  P2 

(2.15) 

+  £Pi  £h  +  ^  +  £?1  ^p2  J 


Since  £  and  £  are  independent  in  themselves  and  £  is  independent  of 
i  h  q 

£p,  eq.  2.15  reduces  to 


y^l2 


^  V 


-  CT"  2 
q 


(2.16) 


The  components  of 


£  and 
*1 


are  derived  in  Appendix  A  as 


£  -  -  Afl,+  A02  (A.  12) 

P1 

and  £  =-AG2-AI^2  (A.  8) 

2 

where  £0,  ,  AG  and  AK  are  angular  changes  defined  in  Appendix  A 
1  2  '4 

(figs.  A-3  and  A-4). 


Therefore , 


S2]  -E[ **i  + 


AG2)(-  4G2  - 
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(2.17) 


or 


€  €  -  E 

L  pi  p2  J  L 


A0.  ^8  +  A®,  A*.,  -  A0, 


A©,  AH. 


The  individual  angular  changes,  A0^>  A 0£  and  A^2»  are  individually 

independent  of  any  other  single  angular  change  so  that  eq.  2.17  may  be 
written  as 


(2.18) 


(2.19) 


Now,  substituting  eqs.  2.11,  2.12  and  2.19  into  eq.  2.13, 


-<  ^22+*r2> 


*  <r2  +  or2) 


(2.20) 


Eqs.  2.11,  2.12  and  2.20  define  all  of  the  parameters  required  in  eq.  2.3. 

In  order  to  continue  the  development  into  Fraser's  method  for  eliminat¬ 
ing  the  dependence  of  the  LOP  errors,  it  is  necessary  to  convert  the  para¬ 
meters  of  eq.  2.3  as  defined  by  eqs.  2.11,  2.12  and  2.19  into  distance  units 
measured  over  a  straight  line  which  is  perpendicular  to  the  exact  and  the 
error  LOP's.  Ultimately,  this  conversion  is  necessary  to  obtain  the  radius 
of  a  particular  circular  probability  density  contour,  which  must  be  in  dis- 

A 

tance  units.  Assume  that  there  is  a  factor  K  which  meets  the  stated  con¬ 
version  requirements.  This  is  expressed  in  equation  form  as 

£,  -  K  €a  (2.21) 

d  a 

where  is  the  error  in  distance  units 

£  is  the  error  in  angular  units. 

8 


and 
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Applying  eq.  2.21  to  £  gives  £  ■  K  £  .  This  implies  that 

l  Id  1  A 

on2-E[(«iV2]  (2-22> 

But.  b[(  SjC/]-  VE[  £!2] 

so  that  (J^  2  =  Kx2  eL2  J  (2-23) 

Similarly,  0^  =  E^  £  2  J  (2.24) 

Using  eq.  2.8  and  substituting  eq.  2.11  or  eq.  2.12  as  appropriate,  eqs. 
2.23  and  2.24  become 

2  -  il2  ( <r2  +  o-2  +  or2 )  (2-25> 

Id  1  i  q  Pj^ 

and  2  =  K22  (  2  +  <772  +  2)  (2-26) 


The  same  procedure  can  be  used  to  convert  the  angular  covariance  term  into 
distance  units.  Using  eqs.  2.14  and  2.21, 

/tl2d  =  E[  <  ^  K2  £2)  J  *  K2  E[  €x  £2  J  (2 .27) 

or,  substituting  eq.  2.19  into  eq.  2.27, 


Note  that  when  the  variances  and  covariance  in  distance  units  given  by 
eqs.  2.25,  2.26  and  2.28  are  substituted  into  the  general  expression  for 
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the  correlation  coefficient 


2.13), 


A 

the  K  factors  cancel  so  that 


//12  _  ^l2d 

<T\  &l  0Td  0~2d 


(2.29) 


The  correlation  coefficient  ^  is  unaffected  by  the  conversion  to  distance 
units  and  is  still  given  by  eq.  2.20.  Chapter  3  is  dedicated  to  the  deri- 

A 

vation  and  analysis  of  the  conversion  factor  K. 

The  variance  and  covariance  terms  in  distance  units,  as  defined  by 
eqs.  2.25,  2.26  and  2.28,  permit  the  application  of  Fraser's  method 
directly.  (Fraser's  method  is  abstracted  herein  for  the  purposes  of  this 
discussion.  The  details  of  the  entire  development  are  given  in  ref.  3.) 
Since  the  LOP  crossing  angle  (  o(  in  fig.  2-1)  is  not  necessarily  90°,  the 
first  step  is  to  transform  the  error  variables  onto  orthogonal  axes.  The 
transformation  is  given  as 


€ 


Id 


and  £  m  Ei  ~  E2  tanol 
2d  sec«C 


(2.30) 


where  E^  and  are  the  error  variables  along  the  new  orthogonal  axes  and 
°Cis  the  LOP  crossing  angle. 

The  transformation  of  eqs.  2.30  can  be  used  to  obtain  the  joint  probability 
density  function  of  E^,  E£  as 
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,  Ei2  2aEiE2 ,  e22  , 

*  _ 1  _ .  +  T  > 


P(  El,  E2)  -  - - — j= 

^o&otV1 

where  or-  «r: 

i  la 


2(1  -  p/>  or/  af;or2  or/ 


(2.31) 


|tal«|  tI  C’7d  '  +  °2d  '  sec‘'°<  -  2/t'l2d  sec°< 


2  2_ 


U12  "  E  [  E1  E2  ]'  ‘  tan*  <  °ld  2  ~/\2d  secoC  > 


/V 


&ld  2  -/\ld  secoC 

^ld  ifo  2  +  ^2d  2  sec  ^  -  $^l2d  sec*. 


Although  the  error  axes  are  now  orthogonal,  the  cross  product  term 
in  the  exponent  of  e  in  eq.  2.31  indicates  the  dependence  of  upon  E^. 

In  order  to  eliminate  the  cross  product  term,  the  coordinate  axes  are  rotated 
by  another  linear  change  of  variables.  The  transformation  is  stated  as 


follows : 


=  x  cos f  -  y  sin./* 


E„  =  x  sin /  +  y  cos/ 


(2.32) 


where  /  is  the  angle  through  which  the  axps  are  rotated. 

The  derivation  of  ref.  3  showr  cnat  the  expression  for  /  must  be 


(2.33) 


2  j°l  CfEl 
&l[  2  ~  2 


This  transformation  for  the  rotation  of  the  coordinate  axes  is  discussed 
in  ref.  4  also.  When  eqs.  2.32  and  2.33  are  substituted  into  eq.  2.31, 
the  resulting  probability  density  function  is 


1  ,  x 


P(x,y) 


2  it  cr^  of 


e 


(— +  £ — ) 

'  0"~  2+  a~  2 

x  y 


(2.34) 


where  or-  Ja^~2  cos  2f  +  072  2  sin2/  +  2U12  sin/  cos/ 


(2.35) 


and 


2  sin2/  +  2  cos2/  -  2U12  sin  f  cos/ 


Eqs.  2.35  define  the  standard  deviations  and  0^  of  the  independent 

x  y 

variables  x  and  y.  Now,  fig.  2-3  can  be  entered  with  the  smaller  of  the 

ratios  or  •  This  value  and  the  desired  accuracy  probability  define 

Oy 

a  value  of  K  .  The  radius  of  the  equivalent  circular  probability  density 
contour  is  calculated  by 

R  -  K  (7£  (2.36) 

where  (?£  is  the  larger  of  and  0 

This  concludes  the  development  of  the  basic  method  for  determining  the 
positioning  accuracy  of  the  method  of  horizontal  sextant  angles.  As  stated 
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TO  CIRCLE 


"IL* 


mm w 
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previously,  the  derivation  and  analysis  of  the  conversion  factor  K  used 
to  convert  the  angular  error  to  distance  error  is  presented  in  Chapter  3. 
The  derivation  of  a  means  for  calculating  the  LOP  crossing  angle  oc  for  a 
three-object,  two-LOP  position  fix  is  included  as  Chapter  4.  Finally,  an 
example  of  a  practical  application  of  the  complete  method  to  a  buoy  posi¬ 
tioning  problem  is  given  in  Chapter  5.  It  should  be  noted  that  the  presen 
tation  will  continue  through  the  example  problem  without  the  benefit  of 
suitable,  published  human  factors  error  data.  The  consequent  assumptions 
made  in  Chapter  5  for  the  purposes  of  the  example  are  explicit  therein. 


fjjjjjfiiiVriiMr  *  •--1  —  - 
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3.  THE  GEOMETRIC  MULTIPLICATION  OF  PRECISION  (GMOP) 

3.1  The  Role  of  the  GMOP 

The.  development  of  Section  2.3  presumed  the  existence  of  a  conversion 

A 

factor  K  ,  defined  by  eq.  2.21,  which  can  change  the  angular  variances  of 
the  individual  LOP's  and  the  angular  covariance  of  the  combined  LOP's  into 
distance  units.  Further,  in  order  to  conform  to  the  geometric  model  of 
fig.  2-1,  it  is  necessary  that,  in  the  neighborhood  of  the  observers'  posi¬ 
tion  P,  the  resulting  distance  error  is  measured  over  a  straight  line  which 
is  perpendicular  to  the  exact  and  the  error  LOP's  (assumption  "d". 

Section  2.1).  It  is  the  purpose  of  this  chapter  to  develop  and  study  a 
suitable  conversion  method  within  the  geometric  constraints  of  the  system. 
The  ultimate  goal  is  to  find  some  function  which  is  independent  of  the  base¬ 
line  length  (the  distance  between  the  fixed  objects,  designateJ  A)  such 
that  the  product  of  this  function,  evaluated  at  any  point  ?,  and  A  will 

A 

yield  K  .  Combining  this  concept  with  eq.  2.21,  this  function,  named  the 
Geometric  Multiplication  of  Precision  (GMOP),  is  expressed  in  equation 
form  as 


GMOP  =  — 
A 


€d 
A  £a 


(3.1) 


where  is  the  error  in  distance  units, 
€  is  the  error  in  angular  units, 
and  A  is  the  baseline  length 


The  reader  familiar  with  electronic  triangulation  systems  will  recognize 
the  analogy  between  the  GMOP  and  the  GDOP  (Geometric  Dilution  of  Precision) 
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commonly  used  for  Che  electronic  systems.  The  significant  difference  be¬ 
tween  the  two  is  that  GMOP  equals  zero  at  the  fixed  objects  0  and  Q. 

At  this  point  a  brief  discussion  of  the  direction  of  the  following 
development  is  appropriate.  The  direction  of  the  development  is  set  by 
the  desire  to  describe  the  GMOP  as  a  function  of  the  magnitude  of  the 
gradient  of  .  The  logic  of  this  choice  is  apparent  when  the 

the  desired  features  of  the  GMOP  are  compared  with  the  inherent  traits  of 
the  gradient  of  .  The  magnitude  of  V/3  converts  distance  traveled 

between  contours  of  constant  ^  (the  LOP's)  along  the  path  of  maximum  change 
in  p  to  an  angular  change  in  £  .*  The  direction  of  along  the  path 

of  maximum  change  in  ^  requires  that  it  is  perpendicular  to  the  contours  of 
constant  .  Fig.  3-1  illustrates  the  relationship  between  the  contours  of 
constant  ^  and  curves  of  V  ■  (The  circular  nature  of  the  curves  of 
as  presented  in  fig.  3-1  is  shown  to  be  accurate  during  the  analysis  of 
Appendix  B.)  Obviously,  the  inverse  of  l^s|  meets  the  requirement  to 
convert  angular  measure  to  distance  units  within  the  geometric  constraints 
of  the  system.  Further,  the  fact  that  y/8  is  perpendicular  to  contours  of 
constant  ^  indicates  that  the  requirement  to  measure  the  error  in  distance 
units  along  a  straight  line  which  is  perpendicular  to  the  exact  and  error 
LOP's  may  be  approximated  by  measuring  the  error  along  curves  of  T  ■ 
Therefore,  there  is  a  strong  possibility  of  fulfilling  all  of  the  require¬ 
ments  of  GMOP  with  a  function  of  |  ^3  |  .  The  development  begins  with  the 
derivation  of  the  mathematical  expression  for  the  LOP  in  Section  3.2. 

3.2  The  Mathematical  Expression  for  an  LOP 

Consider  the  circular  LOP  shown  in  fig.  3-2.  Again,  the  observer's 
position  is  P  and  0  and  Q  are  the  fixed  objects  used  to  generate  the  hori-/ 


*The  conversion  is  from  distance  units  to  radians. 


FIG.  3-2  CIRCULAR  LOP  SHOWING  PERTINENT  PARAMETERS 
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zcntal  sextant  angle.  The  coordinate  system  is  chosen  with  point  0  at 
the  origin.  (It  is  very  important  to  understand  that  this  representa¬ 
tion  and  the  analysis  which  follows  applies  to  any  sextant  angle  generated 
from  two  fixed  objects,,  i.e. ,  the  points  0  and  Q  are  representative  of  any 
set  of  baseline  objects.)  It  is  expedient  in  this  case  to  use  the  polar 
coordinate  system  where  r  is  the  position  vector  of  any  point  in  the  plane 
of  the  coordinate  system.  The  position  vector  r  has  magnitude  r  and  a  direc¬ 
tion  defined  by  the  angle  6  between  r  and  the  baseline  axis.  The  variables 
A  and have  been  defined  as  the  distance  between  0  and  Q  along  the  base¬ 
line  and  the  exact  measurement  of  the  horizontal  sextant  angle  respectively. 
Referring  to  fig.  3-2  and  using  the  law  of  sines ; 


sin  K  sin 


(3.2) 


Substituting  for 


sin  £  fT  -  (  p  +  0)]  sin^g 
Noting  that  sin  £  TP  -  (  ^3  +  ®)J  “  sin  (  f3+  ®) » 


sin  ( ^3  +  0)  sin ^3 

Finally,  sin  (  ^  +  0)  *  sin^cos0  +  cos^i  sin0 
so  that 

r  *  A(cos0  +  cot^gsinO) 


(3.3) 


Eq.  3.3  represents  the  LOP  in  polar  coordinates  in  terms  of  r,  A,  0  and 
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3.3  The  Derivation  of  the  GMOP 

Eq.  3.3  can  be  manipulated  to  yield  an  expression  for  ^  in  terms  of 
r,  A  and  0.  The  manipulation  yields 


*  cot  ^  |^csc0  -  cotoj 

The  general  form  of  the  gradient  of  ^  in  polar  coordinates  is 


(3.4) 


(3.5) 


where  £  and  $  are  unit  vectors  in  the  r  and  9  directions  respectively. 
The  components  of  eq.  3.5  may  be  derived  from  eq.  3.4  as 


cs 


c0sin^3 


and 


-i-  cscQsin2/?  Sln!  ~  cot&-c°sQn 
A  '  cos0  +  cot^sinO 


(3.6) 


(The  details  of  the  derivation  from  eq.  3.4,  through  eq.  3.8  are  presented 
in  Appendix  C.  Also,  an  understanding  of  the  development  through  eq.  3.8 
should  permit  the  reader  to  follow  the  analysis  of  Appendix  B,  the  proof  of 
the  acceptability  of  assumption  "d"  of  Section  2.1.) 


The  magnitude  of  v/s  is  given  by 


(3.7) 


Substituting  eqs.  3.6  into  eq.  3.7  yields 
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m 


_ csc6sin/5 

A(cos6  +  cot^sinfi) 


(3.8) 


Now,  the  GMOP  may  be  derived  from  eq.  3.8  in  the  desired  form  of  eq.  3.1 


sinScsc^  (cos©  -I-  cot^d  sin©) 


(3.9) 


3.4  The  Nature  of  the  GMOP 


At  this  point  the  behavior  of  eq.  3.9  over  the  entire  ranges  of 
and  ©  is  of  interest.  The  equation  may  be  put  in  a  more  suitable  analytical 


form  by  first  expanding  as 

GMOP  =  csc^g  sindcos©  +  csc^rfcot^ sin2© 


(3.10) 


Substituting  into  eq.  3.10  the  trigonometric  identities 
sin© cos 6  =  isin20 


2  11 

and  sin  ©  =  i  -  ■=cos20 
2  2 


yields 


GMOP  =  icsc^sin20  +  icsc^cot^  -  icsc^  cot^ cos2©  (3.11) 


Letting 


C  =  icsc^3  and  D  =  icsc^cot^, 


GMOP  =  D  +  Csin2©  +  (-D)cos20 


or  GMOP  =  D  -  -/C2+D2  cos(2©+y) 


(3.12) 


where  Y  ■  tan  *  (-— ) 
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Let  E 


■/ 


CL  +  D 


and  notice  that  -  =  tan /3  >so  that 
D  ' 


GMOP  -  D  -  Ecos(20  +  p  ) 


(3.13) 


The  nature  of  the  GMOP  is  readily  discernible  from  eq.  3.13.  Particu¬ 
lar  features  include  the  following: 

a.  For  any  £}  ,  GMOP  consists  of  a  constant  D  from  which 
some  portion  of  a  cosine  function  is  subtracted. 

b.  The  sinusoidal  portion  of  GMOP  cycles  once  each  "fT 
radians. 

c.  Increasing  p  between  0°  and  90°  causes  a  decrease 
in  D  and  E  and  shifts  the  cosine  curve  to  the  left. 

d.  Noting  that 


E  «•  -^2  +  j)2  _  ^csc2^  +  icsc^  cot^  «  icsc^  *\J  1  +  cot^3 


l  2  . 
±csc 


eq.  3.13  can  be  written  as 


GMOP  =  icsc^fl  cos^g  -  icsc^3  cos (20  +/*> 


Values  of  0  may  vary  from  gero  for  P  located  at  Q  to  “TT  -  ^ 
for  P  at  0  (see  fig.  3-2)., 

At  these  boundary  values  of  0, 

0  ■  0,  GMOP  »  icsc^?  cos^ -  icsc^  cos ^  -  0 
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6  ■  TT-^  ,  GMOP  «  icsc^  cos p  -  icsc^  cos(2*ff -2^+^)  »  0 


This  is  Che  expected  result  since  the  contours  of  constant  p 
(the  circular  LOP's)  converge  at  0  and  Q  (see  fig.  2-2). 

e.  The  variation  of  0  from  0  to  Tt '  ~/3  causes  the  truncation 

of  the  cosine  curve  at  a  distance  p  from  20  + ^  =  0  and 
from  20  *  2*tP  ;  thus,  all  possible  values  of  20  + 

yield  cos(20  +yfl) ^cos^  so  that  the  GMOP  is  always  positive. 

f.  Intuitively,  the  points  ^  ,  which  lie  on  the  perpendicular 
bisector  of  OQ  in  fig.  3*3,  should  be  the  points  at  which 


GMOP  is  maximum  for  any  given  p  .  As  indicated  in  the  figure, 

•ft  -i 
2 

checked  for  agreement  as  follows: 


for  the  points  ,  0  -  . The  derived  equation  may  be 


^GMOP 


09 


P 


2Esin(20  +  £  )  m  csc^j  sin(20  +/g  )  *  0 


V5 


(3.14) 


for  the  non- trivial  case,  eq.  3.14  implies  that 


20  + 


=  nlT  (n*0,  1,  2,...) 


For  any ^  in  the  range  of  interest,  the  only  one  of  the 
values  which  20  assumes  is  T?  .  Consequently, 

20  +  £  =  IT 

which  yields  0  =  for  maximum  GMOP  ( p  =  const) 

g.  Finally,  the  combined  implications  of  "e"  and  "f"  above  indi¬ 
cate  that  the  values  of  GMOP  are  symmetrical  about  the  perpendicu¬ 


lar  bisector  of  OQ. 
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| 


Figure  3-3.  Circular  LOPs  showing  points  of  maximum  GMOP. 
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Items  "d"  through  "g"  above  indicate  that  the  mathematical  representa 
tion  of  the  GMOP  behaves  as  expected  at  several  check  points.  It  is  con¬ 
cluded  from  these  favorable  results  that  eq.  3.9  does  describe  a  function 
which  satisfies  the  form  requirement  of  eq.  3.1  and  converts  angular  errors 
to  distance  units  within  the  system  constraints.  Fig.  3-4  shows  curves  of 
GMOP  vs  0  for  2  0- degree  increments  of  from  5°  to  125°  (within  the  range 
of  the  newest  sextants).  The  curves  are  presented  to  demonstrate  the 
changing  shape  as  ^  varies.  Additionally,  contours  of  constant  GMOP  are 
shown  in  fig.  3-5.  This  type  of  representation  usually  aids  in  the  visuali 
zation  of  the  physical  situation. 


3.5  The  Feasibility  of  Linear  Interpolation  for  GMOP 

An  investigation  of  the  feasibility  of  linear  interpolation  between 
the  curves  of  fig.  3-4  to  obtain  values  of  GMOP  for  intermediate  values  of 
!$  and  between  curves  of  GMOP  vs  for  incremental  0  in  order  to  obtain 
GMOP  for  intermediate  values  of  0  is  of  interest.  For  the  first  case, 
values  of  GMOP  and  GMOP  (0j,  ^3)  would  be  obtained  from  fig.  3-4 

Then,  the  linear  interpolation  could  be  performed  to  obtain  GMOP  (0^ , ^  2) • 
The  second  case  implies  that  GMOP  (9].>y5i)  and  GMOP  (0-j,  ^ could  be 
used  to  find  GMOP  (02»  1)  by  linear  interpolation.  The  combination  of 

these  two  possibilities  would  yield  a  relatively  concise  table  from  which 
GMOP  could  be  determined  for  any  possible^  and  0.  The  advantages  of 
such  a  presentation  are  obvious.  For;  exact  linear  interpolation  in  the 


first  case 


GMOP 


P 

tionship  between 


fi 


must  be  constant,  i.e.,  there  must  be  a  linear  rela- 
and  GMOP  for  any  constant  0.  Likewise,  exact  linear 


0 


GMOP 


Figure  3-4b.  Curves  of  GMOP  vs.  9.  for  £=  45°,  65°,  85°,  105°,  125°. 
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Figure  3-5.  Contours  of  Constant  GMOP  Drawn  to  Scale 
(Baseline  OQ  of  Length  A) 
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interpolation  in  the  second  case  requires  that 


3gMOP| 


is  constant. 


Te  |  fi 

Therefore,  the  behavior  of  each  of  these  quantities  is  an  indication  of  the 
accuracy  which  can  be  expected  from  linear  interpolation  for  their  respective 


cases.  The  equation  of .  -y 


e 


is  derived  readily  from  eq.  3.9.  The  de¬ 


tails  of  the  derivation  are  shown  in  Appendix  D.  The  result  is 


3GM0P 

30  e 


-Gcsc 


£  sin^cos( £  +  0)  +  2gJ 


(3.15) 


where  G  ■  sinfl 

Note  that  the  quantity  in  brackets  is  a  relatively  "tame"  function  whose 

individual  components  never  exceed  a  value  of  2.0.  The  effect  of  this  quantity 

3 

is  to  modify  the  dominant  -esc  function.  The  degree  of  modification  will 

depend  upon  the  particular  values  of  and  6.  The  coefficient  G  exerts  a 

'N  fZMOPl 

more  commanding  influence  of  the  shape  of  - 1  in  that  the  effect  of 

1 0 

the  constant  0  value  is  not  diluted  by  the  varying  ^6  values.  Curves  of 


3GM0P 


3/3 

nearly  constant) 


0 


are  shown  in  fig.  3-6.  Although  each  curve  exhibits  constant  (or 


3GM0P 


* 


over  some  range  of  ^3  ,  these  ranges  vary  in  size 


and  position  for  different  values  of  0.  In  general,  linear  interpolation 
between  curves  of  GMOP  vs  0  for  constant  incremental  ^  is  acceptable  only 
over  limited  ranges  of  ^  which  differ  with  changing  0. 


The  equation  for 


3  GMOP 


3  0 


is  taken  from  eq.  3.14  as 


3  GMOP 


3e 


=  2Esin(20  +^3  ) 


(3.16) 


,  2 

where  E  =  esc  , 
2 
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Figure  3-7. 
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For  any  constant^,  eq.  3.16  represents  a  sine  curve  of  amplitude  csc^S. 
The  most  acceptable  ranges  of  0  for  linear  interpolation  purposes  are  those 
which  include  the  section  of  the  curve  near  a  peak.  Since  the  phase  angle 
of  the  curves  changes -with  0  (the  phase  angle  is,  in  fact,  equal  to ^  ), 
these  acceptable  ranges  will  vary  for  varying  ^3  .  Additionally,  the  ampli¬ 
tude  of  the  curves  varies  inversely  with^J*j  and,  since  a  smaller  amplitude 
will  effect  a  "flatter"  curve,  the  range  of  0  over  which  linear  interpola¬ 
tion  is  acceptable  increases  with  increasing  £0.  Curves  of  are 

presented  in  fig.  3-7.  The  results  of  figs.  3-6  and  3-7  prompt  a  general 
caution  against  linear  interpolation.  Until  the  acceptable  error  in  GMOP 
is  determined  (as  the  error  effects  the  calculation  of  ^"ld,  ^2d  and//<l2d) 
it  is  advisable  to  calculate  the  required  value  for  any  specific  case 
directly  from  eq.  3.9. 


*for 


OJy as 


90c 


-63- 


4.  THE  LOP  CROSSING  ANGLE  («<) 

4.1  The  Calculation  of  oC 

This  chapter  is  dedicated  to  the  derivation  and  presentation  of  a 
method  for  calculating  the  LOP  crossing  angle  cC  as  an  input  to  the  basic 
method  of  determining  the  positioning  accuracy  developed  in  Section  2.3. 
The  LOP  crossing  angle  is  defined  as  the  acute  angle  between  the  tangents 
to  the  LOP's  at  the  observer's  position  P  (see  fig.  2-1).  It  is  shown 
physically  for  the  respective  cases  of  figs.  1-4  in  figs.  4-1.  The  effect 
of  o(  on  the  reliability  of  a  position  fix  will  not  be  addressed  speci¬ 
fically  in  this  report.  (The  deterioration  of  fix  reliability  with  the 
deviation  of  from  90°  is  demonstrated  effectively  in  fig.  7  of  ref.  2.) 

Consider  the  illustration  of  one  LOP  in  fig.  4-2.  The  development 
of  Section  1.1  provides  the  value  of  ^COQ  ■  (point  C  is  the 

center  of  the  circular  LOP).  It  follows  from  fig.  4-2  that 


^  -  8  -  £C0Q  =  0  -  -LL_  +  ^  (4.1) 

Since  ^OPC  is  isosceles  (two  of  its  sides  are  radii  of  the  circle), 
<£P0C  =  a»0PC;  and  in  fig.  4-2  is  given  by 


(4.2) 


This  expression  for  ijj  can  be  derived  readily  for  configurations  encompas¬ 
sing  all  possible  combinations  of  jS  and  0. 

The  two  general  configurations  of  the  tangents  to  the  LOP's  and 
the  position  vectors  "x  at  the  observers'  position  P  are  shown  for  the  case 
of  fig.  4-la  in  figs.  4-3.  From  fig.  4-3a, 


TANGENT  TO  LOP  AT  P 


Figure  4.2  One  LOP  showing  the  angle  ♦ 
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«t-  'h  '<fi  -/3i> 


(4.3) 


Substituting  eq.  4.2  into  eq.  4.3  yields 


°C  =  02  +  p  2  '  ®1 

In  like  manner,  from  fig.  4-3b, 

°C *  •n,-['p2  -  <j*i 

or  =  *TT  -  (02  +/32  -  0X) 


(4.4) 


(4.5) 


Figs.  4-4  show  the  two  general  configurations  for  the  case  of  fig.  4-lb. 
From  fig.  4-4a, 


«■-  fi  -  Vi> 

or  -  0i  -  G2  -  ^2  “  "  (02  +  ^  2  _  el) 

From  fig.  4-4b, 

AC  -  TT-  |^1  -  ( f*2  +/3i)J 

or  =  'TT  -  (©i  -  @2  "  P  2^  “  'K'_  (®2  +  ^  2  "  ®l^j 


(4.6) 


(4.7) 


Eqs.  4.4,  4.5,  4.6  and  4.7  can  be  used  to  form  the  general  rule  for 
determining  the  LOP  crossing  angle  c£,  stated  as  follows: 


=  02  +  ^2  ”  ®1  >|  ®2  ”  ®1  |  ^  90® 

o(_  *  180°  -  ©2  +  ^  2  “  ®1  I >  2  ”  ®1  ^  90° 


(4.8) 


where  the  absolute  value  of  (0£  +  ^2  "  ®l)  ta^en  to  allow  for  the 
negative  values  of  this  quantity  which  are  encountered  for  the  case  of 
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fig  4-lb  (eqs.  4.6  and  4.7). 

It  should  be  noted  that  when  ■  0,  the  two  horizontal  sextant  angles 
define  the  same  LOP  (fig.  4-5).  This  is  called  commonly  a  swing  angle 
and  is  characterized  by  02  +  ^2  =  ®1»  The  disastrous  implications  of 
selecting  objects  which  create  this  situation  are  apparent  from  fig.  4-5. 
The  proper  sextant  angles  can  be  measured  from  any  point  on  the  LOP. 

Also,  note  that  eqs.  4.8  express  0(  in  terms  of  @2>^2  and  9],.  Since 
these  values  are  required  for  the  calculation  of  GMOP  (eq.  3.9),  the 
determination  of  involves  very  little  additional  effort. 

The  developments  of  this  chapter  along  with  those  of  Chapters  2  and 
3  form  the  complete  method  for . evaluating  the  accuracy  of  a  position 
fix  obtained  from  horizontal  sextant  angles.  Since  these  developments 
are  cluttered  somewhat  by  the  supporting  theory,  an  example  problem  is 
provided  in  Chapter  5  to  demonstrate  the  application  of  the  complete 
method  to  a  practical  situation. 

Additional  Note:  In  order  to  apply  rigorously  the  method  for  deter¬ 
mining  positioning  accuracy  developed  herein,  the  proper  absolute 
value  and  sign  of  the  crossing  angle  o(  must  be  used.  The  absolute 
value  of  is  determined  from  eqs.  4.8.  If  *  is  measured  from  LOP-2 
to  LOP-1  in  a  counter jlockvise  direction  (figs.  4-3b  and  4-4a) ,  the 
sign  is  positive;  if  measured  from  LOP-2  to  LOP-1  in  a  clockwise 
direction  (figs.  4-3a  and  4-4b),  the  sign  is  negative.  However,  the 
sign  of  *  affects  both  Ui^  (as  defined  for  eqs.  2.31)  and  J*  (eq.  2.33) 
directly.  The  manipulations  of  eqs.  2.35  neutralize  the  sign  effects, 
ie.,  the  values  of  ^  and  CJ.'  calculated  from  eqs.  2.35  are  unaffected 
by  the  sign  of*..  Consequently,  to  simplify  the  mechanics  of  applying 
the  method,  only  positive  s  are  used. 
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5.  AN  EXAMPLE  PROBLEM 
5.1  The  General  Procedure 

The  developments  of  the  preceding  chapters  can  be  combined  now  with 
the  concepts  of  Appendix  A  to  formulate  the  general  step-by-step  procedure 
for  determining  the  accuracy  of  a  horizontal  sextant  angle  position  fix. 
The  steps  of  the  procedure  are  listed  as  follows: 

Step  1:  Determine  the  angular  variances  of  the  individual  LOP's  and  the 
angular  covariance  of  the  combined  LOP's  using  eqs.  2.11,  2.12  and  2.19. 


CT?- 

or2  +  ct'2  +  er' 2 

l  9  Pj 

(2.11) 

or2  +  cr2  +  err2 

h  w  q  p 

2 

(2.12) 

^12  “ 

-  <rV  *  < 

(2.19) 

The  definition  of  the  component  variances  can  be  derived  from  the  relation¬ 
ships  established  in  chapters  1  and  2  and  Appendix  A.  Consider  first  the 
variance  .  The  angular  error  is  given  in  Chapter  2  as 

-  ec  +  er  <5.» 

As  defined  in  Section  1.3,  £0  consists  of  the  error  in  placing  the  arm  PO 
directly  through  the  center  of  point  0  (fig.  1-5)  and  the  error  in  charting 
the  position  of  the  fixed  object  0.  For  the  purposes  of  this  example,  it 
is  assumed  hereafter  that  the  charted  position  of  each  of  the  points  0,  Q 
and  Q'  is  exact.  Therefore,  £,Q  consists  solely  of  the  error  in  placing 
the  arm  PO.  The  error  Gr  is  defined  in  Section  1.3  as  the  combination 


of  the  independent  system  errors  including  the  three-arm  protractor  instrument 
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error  and  the  sextant  instrument  and  human  errors.  In  equation  form, 


€  =  £,  +  C  +  £ 

r  TJa  sti  sth 


(5.2) 


where  £ is  the  three-arm  protractor  instrument  error;  Cgti  is  the 
sextant  instrument  error;  and  £sth  is  the  sextant  human  error. 

Eq.  5.2  can  be  substituted  into  eq.  5.1  to  give  £^  as 


£ i  *  £„  +  £ja  +  £sti  +  £s 


(5.3) 


Since  each  of  the  components  G0,  £33,  £stj.  and  £gth  are  individually 

independent,  eq.  5.3  can  be  used  in  eq.  2.8  to  yield 


of  -  e[  e,2]  -  <r?  +  <7?a2  +  oTti2  +  oTn, 


(5.4) 


Similarly, 


cr2  -  e|"  e, 2 1  =  crt2  +  <jr2  +  cr; 2  +  cr  2  (5.5) 

v  h  L  h  J  ‘l  3a  sti  sth 

9  _ *  2 

The  variances  and  can  be  determined  by  using  the  expres- 

P1  P2 

s^jns  for  £  and  £  given  in  Appendix  A  as 


€p  -  -  A01 


£  --A92  -  M2 


(A. 12) 


(A. 8) 


in  eq.  2.8.  Since  each  of  the  angular  changes  A©2  anc*  ^^2  are 

independent  of  any  other  single  angular  change,  the  results  ar» 


cr2  =  e  I”  e  2 1  =  or 2  +  2 

L  Pi  J  ^2  A91 


(5.6) 


and  CT'  2  =  E  f  £  2  =  <ZT 

P2  L  P2  J  A02 


2+  ctT'2 


(5.7) 
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Now,  substituting  eqs.  5.4  through  5.7  into  eqs.  2.11  a-.d  2.12, 


q2  -  <r02  +  oj,2  +  <^!2  +  xJ  +  cf  +  err  2  +  cqs 


(5.8) 


cr22  ■  o;, 2  +  oj  +  0Ttl2  +  <q:h2  +  crq2  +  ^  +  o^2  (5.9) 


Eqs.  5.8,  5.9  and  2.19  show  completely  the  component  variances  of  » 

2 

(T^  and//f'2  which  will  be  considered  for  this  example  problem.  The  calcu¬ 
lation  of  the  values  of  the  component  variances  is  considered  in  Section 


A  A  _ _ _ _ 

Step  2;  Determine  and  1^,  the  factors  required  to  convert  O^,  CT^ 
and^t^2  from  angular  units  (radians)  to  distance  units.  Combining  eqs. 


3.1  and  3.9, 


£  =  (A)GMOP  =  AsinOcsc^  (cos9  +  cot ^  sin8)  (5.10) 

where  A,  9  and  j$  are  shown  in  fig.  3-2. 

If  the  values  of  A^,  9^  and  ^5^  for  LOP-1  are  substituted  into  eq.  5.10, 

A  A 

the  value  of  is  obtained;  using  A^,  0£  and 2  for  L0P-2  produces  l^* 

A 

It  should  be  noted  that  hand  calculations  of  K  can  be  facilitated  by 


recognizing  that 


r  =  A(cos0  +  cot^ysin9) 


so  that 


K  =  rsin9csc^ 


(3.3) 


(5.11) 


where  r  is  the  distance  from  0  to  P  for  LOP-1  and  from  Q  to  P  for  L0P-2. 


The  variances  and  the  covariances  are  converted  to  distance  units  as 


follows: 
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crr '  =  k  cr 

Id  11 


K2  ^2 


(5.12) 


/^12d  K1  K2  /*12 


Step  3;  Calculate  the  LOP  crossing  angle  from  eqs.  4.8 

|e2+/32  -«l|*  IVA  -  ®1  I1  900 
«-  180°  -U  +/a2  -  .  hv2-ei|>90° 


(4.8) 


Step  4:  Use  «  C7^j»  y^2d  from  eqs‘  ^,12  and  from  eqs*  /*-8  in  the 

expressions  for  cC  ,  OT^  ,  U  and  /?  as  they  are  defined  for  eq.  2.31. 

^2  12  j  1 


\'<ru 


%  '  V^2  +  ^  '  2/t12d  8eC°t 


^12  *  ”  tanct  <0Td2  '/^2d  SeCot  > 


(5.13) 


°Td  ‘/fed sec* 


^  °*1  S  CTTdV  +  ^  •  2/<12d 


seco<_ 


Step  5:  Calculate  the  angle  of  axis  rotation  from  eq.  2.13. 


/.  f  tan-1 


2/*l  ^1  ^2 

or2  -  cr2 


(2.33) 
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Step  6:  Calculate  CT'  and  cT'  by  using  the  values  of  O'"  .  CT"  and  U 

x  y  EX  e2  12 

obtained  from  eqs.  5.13  and  the  value  of  f  from  eq.  2.33  in  eqs.  2.35. 


(2.35) 


Step  7:  Enter  the  curves  of  fig.  2-3  with  the  smaller  of  the  ratios 

CTv .  .  _  _ 


or  and  the  desired  accuracy  probability  to  obtain  a  value  of  K.  The 
radius  of  the  circular  probability  density  contour,  whose  center  is  the 
center  of  the  black  dot  used  to  mark  the  charted  buoy  anchor  position, 
within  which  the  particular  horizontal  sextant  angle  position  fix  can  be 
expected  to  be  for  the  particular  accuracy  probability  is  calculated 


using  eq.  2.36. 


R  =  K  Of  (2.36) 

b 

where  (JT  is  the  larger  of  (7"  and  £7'". 
b  x  y 

It  can  be  seen  that  hand  calculations  of  the  positioning  accuracy 
for  any  particular  problem  are  relatively  long  and  tedious.  However,  the 
method  does  lend  itself  nicely  to  computer  programming,  where  the  para¬ 
meters  peculiar  to  the  problem  would  be  the  required  inputs.  Many  calcu¬ 
lations  could  be  handled  by  a  computer  in  a  short  time.  Nevertheless,  a 
sample  hand  calculation  will  be  carried  out  in  Section  5.3  to  illustrate 
the  mechanics  of  using  the  step-by-step  procedure. 


-76- 


5.2  Probabilistic  Variable  Definitions 

This  section  addresses  the  problem  of  defining  the  means  of  calculat¬ 
ing  the  component  variances  of  C7^2.  U andy^l2’  shown  in  eqs*  5'8’ 

5.9  and  2.19.  The  reader  is  referred  to  Appendix  A  for  the  derivation  of 

the  equations  which  express  €q.  6q,»  A  >  A62  and  in 

terms  of  measurable  distance  errors  in  placing  the  three-arm  protractor 
arms  directly  through  the  centers  of  0,  Q  and  Q'  and  in  placing  the  pivot 
point  directly  over  P.  Applying  eqs.  2.8  and  the  implication  of  eqs.  2.23 
and  2.24  to  eqs.  A.l,  A. 2,  A.3,  A. 6,  A.9  and  A. 10  yields 


(5.14) 


2 

2  _  CTs/  , 


2 
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2 

cm 

2 
r 

3 

where  r  ,  and  are  the  lengths  PO,  PQ  and  PQ*  respectively.  Limited 
experimental  results  from  twelve  samples  of  each  variable  taken  by  two 
subjects  indicate  that  it  is  reasonable  to  assume  the  following: 

a.  The  distributions  of  y^,  y  and  y^  are  identical.  Each  approxi¬ 
mates  a  uniform  distribution  which  is  symmetrical  about  the 
center  of  the  black  dot  used  to  chart  the  position  of  the 
particular  fixed  object.  Further,  it  seems  that  the  human 

mind  can  edit  possible  error  positions  of  0,  Q  and/or  Q'  which 
lie  beyond  the  limits  of  a  circle  whose  center  is  that  of  the 
black  dot  but  whose  diameter  is  equivalent  to  that  of  the 
smallest  black  dot  used  to  chart  the  positions  of  fixed  objects 
ashore  (approximately  0.01  inch  on  all  NOS  charts). 

b.  The  error  points  P£  (figs.  A-3  and  A-4)  have  a  circular 
uniform  distribution  about  the  center  of  the  black  dot  used 
to  chart  the  buoy  anchor  position.  Again,  error  points  out¬ 
side  a  circle  of  diameter  equal  to  that  of  the  small  black 
dot  used  to  chart  the  positions  of  fixed  objects  ashore  seem 
to  be  eliminated.  Consequently,  the  diameter  of  the  circular 
distribution  is  approximately  equal  to  that  of  the  small  black 
dot. 

Therefore,  the  variances  on  the  right-hand  sides  of  eqs.  5.14  are  equal. 

The  standard  deviation  of  the  uniformly  distributed  variables  is  defined  by 

d 

yir 


cr= 


(5.15) 
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where  d  is  the  diameter  of  the  small  black  dot. 

Using  0.01  inch  as  the  value  of  d  and  generalizing  eq.  5.15  for  any  chart 
scale  factor  N, 

CTl  — (ML)N  =  (2.4  x  10"4)N  (5.16) 

(12)  (  Vl2) 

where  CT"  is  in  feet. 

This  expression  for  o'  can  be  used  in  each  of  the  eqs.  5.14  to  determine 
the  indicated  component  variances.  However,  it  must  be  emphasized  that 
eq.  5.16  is  an  expression  for  ^"which  is  based  upon  assumptions  made  from 
very  limited  data.  Its  accuracy  will  remain  suspect  until  it  is  substantia¬ 
ted  by  much  more  experimental  data. 

Preliminary  experimental  results  presented  in  ref.  7  indicate  that 
the  standard  deviation  of  the  human  error  in  measuring  angles  using  the 
sextant  is  approximately  6 '05.5".  It  will  be  assumed  for  the  purposes 
of  this  example  problem  that 

<T"  =  6'05. 5"  =  17.7  x  10-4  radians  (5.17) 

sth 

Assuming  that  the  sextants  used  are  adjusted  perfectly,  the  instrument 
error  in  measuring  angles  with  the  sextant  is  due  to  the  limited  precision 
of  the  instrument.  A  general  assumption  which  is  made  to  approximate  the 
expected  error  due  to  this  limited  precision,  called  the  granularity,  is 
that  the  error  is  uniformly  distributed  between  the  smallest  graduations 
on  the  instrument.  Usually,  sextants  offer  marked  precision  to  the  nearest 
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tenth  of  a  minute.  Therefore, 


cr"  =  .  _ 

Usti  (60)(57.3)(-/l2) 


0.08A  x  10  radians 


(5.18) 


Similarly,  since  the  most  commonly  used  three-arm  protractors  are  marked 
to  two  minutes  of  arc, 


~ — 7=T  =  1.7  x  10"4  radians 


3a  (60)(57.3)(  VlJ) 


(5.19) 


Eqs.  5.1A,  5.16,  5.17,  5.18  and  5.19  define  all  of  the  component  variances 
2  2 

of  (T^  >  CT 2  andy!4^2  f°r  the  purposes  of  the  example  problem. 


5.3  The  Actual  Problem 

The  problem  which  will  be  presented  to  illustrate  the  application 
of  the  method  for  determining  positioning  accuracy  which  is  derived  in  this 
report  is  that  of  positioning  the  anchor  of  Bridgeport  Entrance  Lighted  Bell 
Buoy  // 10.  A  portion  of  C&GS  (now  NOS)  chart  i?220  which  includes  the  buoy 
anchor  position,  labeled  P,  and  the  labeled  fixed  objects  0,  Q  and  Q'  is 
shown  in  fig.  5-1.  The  following  parameters  are  taken  from  the  chart: 


N  =  20,000 


Pi- 


=  92  12' 


=  52  32* 


6,  =  50°06 ' 


02  =  7A  AA ' 


r  =  6A50  ft 
r2  =  8A00  ft 
r3  =  10,050  ft 
A  =  10,800  ft 


8A00  ft 
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Step  1:  Using  eq.  5.16, 

CT-  (2.4  x  10_4)(N)  =  4.8  ft 

Substituting  this  value  into  eqs  5. 14  yields 


CT'2  =  O'?*  2  =  (4‘-  )  ■  55.4  x  10”®  radians 
o  6450 

of  -  C7^22  -  (£t§0>2  -  32.7  x  10-8  radians 

Cr;2  -  <^22  =  <i0/650)  *  22  • 8  x  10"8  radians 

Using  these  values  and  those  of  eqs.  5.17,  5.18  and  5.19  in  eqs.  5.8,  5.9 
and  2.19  yields 


-4 

CT”  *  22.2  x  10  radians 
1 

_ _  -4 

-  20.7  x  10  radians 


-8 

s  "  ^.4  x  10  radians 


Step  2;  Using  eq.  5.10, 

*  5030.5  ft/radian 

^  *  10,237.0  ft/radian 
Eqs.  5.12  can  be  used  to  obtain 

err'..  =  11.2  ft 

Id 

<J~U  *  21.2  ft 

y“i2d  *  -  33.7  ft2 


-82- 


Step  3:  Using  eqs.  A. 8, 
oC  =  76.2° 

Step  4:  Substituting  and  into  e<JS‘  5,13  yields 


CT  =  11.2  ft 
Ei 

O'  =  22.4  ft 
E2 

2 

U  =  -  65.5  ft 


=  -  0.26 

Step  5:  Use  eq.  2.33  to  find 
/-  9.6° 


Step  6:  Calculate  <7-"  and  O' '  using  eqs.  2.35  as 

cr^  =  io.7  ft 


O'"-  22.6  ft 

Step  7:  -£*<J£JL  so  fig.  2-3  is  entered  with 
^  CTy  CT c 

desired  accuracy  probability  of  95%  to  obtain 


=  0.47  and  the 

CT'y 


K  =  2.025 

Using  eq.  2.36,  the  radius  from  the  center  of  the  black  dot  within  which  the 
sinker  can  be  expected  to  be  placed  for  a  95%  accuracy  probability  is 

R  =  45.8  ft 


This  step-by-step  procedure  can  be  used  in  like  manner  to  find  the  position¬ 
ing  accuracy  for  any  other  problem. 
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APPENDIX  A 


The  discussion  of  error  sources  in  Section  1.3  presents  the  errors 
in  measuring  the  horizontal  sextant  angler  using  the  three-arm  protractor 
as  angular  errors  which  modify  the  exact  angular  measurements.  These 
errors  result  from  the  improper  placement  of  the  protractor  arms  PO,  PQ 
and  PQ'  (fig.  1-5)  and  the  improper  placement  of  the  pivot  point  over  the 
point  (observer's  position)  P^  This  appendi:  addresses  the  problem  of  de¬ 
fining  the  statistical  distribution  relationships  of  the  angular  errors 
to  measureable  distance  errors  due  to  the  use  of  the  three-arm  protractor. 
These  distributions  are  necessary  to  determine  some  of  the  variances  and 
the  covariances  used  in  the  development  of  Section  2.3. 

Consider  first  the  angular  errors  CQ ,  and  which  result  from 

the  improper  placement  of  arms  PO,  PQ  and  PQ* .  If  the  protractor  pivot 
point  is  fixed  directly  over  the  point  P  and  the  center  arm  PQ  is  positioned 
through  point  Q,  all  of  the  angular  error  will  be  contributed  by  the  in¬ 
correct  placement  of  the  arm.  Fig.  A-l  shows  the  exact  placement  of  the 
arm  PQ  and  a  hypothetical  incorrect  placement  PQfc  .  The  vector  y2  is 
constructed  from  PQ  to  PQt  ,  perpendicular  to  FQ  at  Q.  Assuming  that  the 
error  point  Qe  is  always  within  the  small  blac’c  dot  with  which  the  posi¬ 
tion  of  an  object  is  charted,  the  magnitude  y£  is  always  much  smaller  than 
the  length  PQ.  Consequently,  £  is  small  and 


e 


q 


tan-1^) 

r2 


21 

r2 


(A.  1) 


where  y^  and  r^  are  shown  in  fig.  A-l  and  £.  is  in  radians. 


*It  is  assumed  that  the  charted  position  of  each  of  the  points  0,  Q  and 
Q*  is  exact. 
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P 


FIG.  A-1  THE  ERROR  IN  PLACING  ARM  PQ 
THROUGH  THE  CENTER  OF  Q 


p«  P 


FIG.  A-2  THE  ERRORS  AFFECTING  DUE  TO  THE 
INCORRECT  PLACEMENT  OF  PIVOT  POINT 
OVER  P 
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Eq.  A.l  implies  that  the  distribution  of  £  can  be  specified  completely 
by  the  one -dimensional  distribution  of  divided  by  Tj.  Note  that  ^ 
is  constant  for  any  particular  problem.  Similarly,  the  distributions  of 
and  '  are  defined  by 


(A. 2) 


(A. 3) 


where  y^  and  y^  are  the  magnitudes  of  the  y"  vectors  from  PO  and  PQ'  to 
PO^  and  PQ^  ,  perpendicular  at  0  and  Q'  respectivelyjand  r^  and  r^  are 
the  lengths  PO  and  PQ1 . 

The  determination  of  practical  expressions  for  £.  and  £,  is  a 

P1  p2 

more  complicated  problem.  If  the  arms  PO,  PQ  and  PQ1  are  fixed  through 
the  centers  of  0,  Q  and  Q' ,  all  of  the  angular  error  is  due  to  the  improper 
placement  of  the  pivot  point  over  point  P.  Fig.  A-2  shows  the  effect  of 
the  error  on  the  measurement  of  ^OPQ.  The  incorrect  placement  of  the 
pivot  point  changes  the  orientations  of  arms  PO  and  PQ  which  causes  the 
angular  changes  AQj  and  Since  the  sum  of  the  three  angles  must 

remain  180°,  the  measured  value  of  the  horizontal  sextant  angle  may  be 


affected  also.  The  angular  effect  on  A  is  designated  £.  ;  in  fig.  A-2, 

/  i  pi 

A 1C  =  A  +  £,  ,  whe~e  the  sign  of  £  depends  on  the  situation. 

rifc  r  i  P]_  Pi 


j3 1  +  €p^  +  0i  +  49i  +  *Vi  +  A/[i 


(A. 4) 


where  the  signs  of  and  at.  also  depend  on  the  situation. 
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But  yd  i  + 


Therefore,  eq.  A. 4  reduces  to 


£pi+  A61+  Afc-0 


6p  =  -(  A©!  +  A^) 


(A.  5) 


Eq.  A. 5  implies  that  the  distribution  of  £^  can  be  derived  from  the 
distributions  of  and  A^  •  The  distributions  of  and 

can  be  determined  using  the  method  by  which  £Q,  and  £^i  were  defined. 
In  fig.  A-3  vector  s*  is  constructed  from  PO  to  P£  0,  perpendicular  to  PO  at 
P.  Since  is  small, 


(A. 6) 


where  r^  is  the  length  PO. 

Similarly,  if  "t  is  constructed  from  PQ  to  P£  Q,  perpendicular  to  PQ  at  P 
(fig.  A-3), 


2 


(A. 7) 


where  ^  is  the  length  PQ. 

The  distribution  for  the  angular  error  in  the  measurement  of  A  QPQ'  due 
to  the  improper  placement  of  the  protractor  pivot  point  can  be  derived  from 
fig.  A-4  as 


£p2  =  AQ2  + 


(A. 8) 
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where  A0„  *  — 

2  r2 


(A. 9) 


and  (A. 10) 

*2  r3 

Notice  that  Z.0QQ'  =  +  Gj.  Since  the  arms  are  fixed  through  the  centers 

of  their  respective  points,  Z.0QQ'  is  constant  and  any  change  in  creates 
an  equal  and  opposite  change  in  Ki  or 

02  -  -  ^  (A. 11) 

Eqs.  A. 7  and  A. 9  reflect  this  relationship.  Substituting  eq.  A. 11  into 
eq.  A. 5, 


(A. 12) 


Eqs.  A. 12  and  A. 8  are  expressions  which  define  the  distributions 

of  £  and  £  .  The  distributions  of  the  components  are  defined  by  eqs. 

P1  p2 

A. 6,  A. 9  and  A. 10.  The  important  implication  of  this  appendix  is  that,  if 
the  one-dimensional  distributions  of  y  ,  y^  and  y^are  symmetrical  about  the 
centers  of  the  black  dots  used  to  chart  the  positions  of  the  fixed  objects 
and  if  the  distribution  of  the  error  points  (figs.  A-3  and  A-4)  about 
the  center  of  the  black  dot  used  to  chart  the  position  of  the  buoy  anchor 


is  circular,  suitable  samples  of  y  ,  y3»  y^,  s1',  t  and  u  taken  from  one 
experiment  can  be  used  to  determine  the  distributions  (and  the  values  of  the 


probabilistic  variables)  of  the  angular  errors  €-q,  and  and  the 

angular  changes  A02  and  &K„,  which  define  £  and  £  .  The 

U  P1  P2 

probabilistic  variables  describing  these  errors  are  necessary  for  the  imple¬ 


mentation  of  the  method  for  determining  positioning  accuracy. 
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APPENDIX  B 

Assumption  "d"  of  Section  2.1,  made  with  reference  to  fig.  2-1, 
is  restated  for  convenience  as  follows: 
in  the  neighborhood  of  the  observers'  position  P, 

d.  the  LOP's  approximate  straight  lines  and  the  angular 
errors  £,  can  be  transformed  into  distance  errors 
measured  along  a  straight  line  which  is  perpendicular 
to  the  exact  LOP  and  the  error  LOP. 

It  is  the  purpose  of  this  appendix  to  present  the  analysis  which  proves 
the  general  acceptability  of  this  assumption.  As  stated  in  ref.  2,  the 
LOP's  approximate  straight  lines  in  the  neighborhood  of  P  if  the  radius 
of  curvature  of  the  exact  LOP  is  much  larger  than  the  one  -  or'  error  (in 
distance  units).  This  requirement  is  represented  in  equation  form  as 


JL»  1  (B-1) 

°d 


where  R  is  the  radius  of  the  exact  circular  LOP  and  is  the  standard 
deviation  of  the  random  error  in  distance  units. 

Fig.  B-l  i.^  a  reconstruction  of  fig.  l-3a  which  includes  additionally 
some  angular  and  linear  measurements  derived  in  Section  1.1.  From  this 
figure  it  can  be  seen  that 


A  1 

R  =  ~y~  A  esc 


(B.2) 


The  combination  of  eqs.  3.1  and  3.9  yields 
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FIG.  B-1  RECONSTRUCTION  OF  FIG.  1-3a  SHOWING 
ADDITIONAL  ANGULAR  AND  LINEAR 
MEASUREMENTS 


£  ( 6 =90°) 

Figure  B-2.  Graphic  Representation  of  Equation  B.5  for  0<  $  <5°. 
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CT^  =  AsinQcsc^  (cos9  +  cot^sin9)C7^ 
Now,  eq.  B.2  can  be  divided  by  B.3  to  give 


(B.3) 


A 

R  _ 


2  sin0(cos9  +  cot^  sin9)  £T^ 


(B  4) 


Preliminary  indications  are  that  <K~  10  minutes;  therefore,  o'- 

3 

-3 

3  x  10  radians  will  be  assumed  for  the  purposes  of  this  analysis.  The 

A 

D 

smallest  values  of  -  will  be  encountered  when  the  denominator  of 

eq.  B.4  is  very  large.  It  can  be  seen  that  the  trigonometric  functions 

of  9  in  eq.  B.4  are  very  "tame",  i.e.,  none  of  them  will  ever  exceed  the 

value  of  1.0.  Since  err  is  considered  constant,  the  volatile  term  is 

a 

cot^d  which  will  approach  infinity  as  ^  approaches  zero.  For  small 
values  of  £5  (0  -  ^  ^  5°),  setting  9  to  90°  will  permit  a  clear  view  of 
the  worst  situations  which  can  be  expected,  identifying  different  areas 
over  the  range  of  ^  from  zero  to  five  degrees  for  which  eq.  B.l  may 
not  be  satisfied.  If  9  =  90°,  eq.  B.4  may  be  rewritten  as 


(cT^)  0  .  n  -  U.7  x  l°2)ta ofi 


(B .  5) 


Fig.  B-2  is  the  graphic  representation  of  eq.  B.5  for  0  -  A  -  5°.  It 

a  ' 

R  > 

can  be  seen  from  this  figure  that,  if  -rp  &  10  is  considered  sufficient 

°d 

to  satisfy  the  implication  of  eq.  B.l,  the  first  part  of  assumption  "d" 
is  acceptable  for  ^  -  3°22' .  It  should  be  noted  that  this  result  implies 
the  acceptability  of  the  assumption  of  linear  LOP's  in  the  neighborhood 
of  P  for  all  possible  values  of  9  corresponding  to  any  ^  greater  than  or 
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equal  to  3°22'  (the  upper  limit  on  is  the  maximum  angle  which  can 
be  measured  by  the  sextant  -  126°  for  the  newest  sextants).  However, 
it  does  not  preclude  the  acceptability  of  this  assumption  for  some 
possible  values  of  0  corresponding  to  ^$'s  less  than  3°22'.  In  fact, 


there  is  some  0  which  yields  ■  —  —  ^  10  for  every  possible  value  of 
except  P  -  0.  The  logical  conclusion  is  that,  in  general,  the  first 
part  of  assumption  "d"  is  acceptable. 

The  second  part  of  assumption  "d"  is  acceptable  if  the  radius 
of  curvature  of  the  particular  curve  of  V  along  which  the  error  is 
measured  is  much  larger  than  the  one  -Of  error  (in  distance  units). 

The  logic  of  this  statement  can  be  seen  intuitively  from  the  representa- 
tion  of  the  contours  of  constant  p  and  the  curves  of  shown  in 

fig.  3-1.  Unfortunately,  the  definition  of  the  radius  of  curvature  of 
the  curves  is  not  readily  apparent  from  the  figure.  In  order  to 

verify  the  shape  and  determine  the  expression  for  the  radius  of  curvature 
of  these  curves,  it  is  most  useful  to  develop  a  scalar  expression,  similar 
to  eq.  3.4,  which  describes  the  whole  family  of  V/3  curves.  To  accom¬ 
plish  this,  the  property  of  the  mutual  orthogonality  of  the  contours  of 
constant  ^  and  the  curves  is  exploited.  By  the  definition  of 

the  gradient,  there  is  some  function  (j>  which  is  constant  on  the  V 
curves  and  for  which  the  contours  of  constant^  are  the  gradient  curves. 
This  concept  is  analogous  to  the  velocity  potential-stream  function  rela¬ 
tionship  of  fluid  dynamics  or  the  potential-flux  relationship  of  electro¬ 
magnetic  theory.  The  fact  that  V^>  •  =  0  everywhere  is  used  in 

the  form  of  the  Cauchy-Rietnann  equations  (ref.  5)  in  polar  coordinates. 


-94- 


These  are  stated  as  follows: 


o)r  r  7© 


and  1  . 

“7®“" 


(B.6) 


Eqs.  B.6  can  be  used  to  find  the  function  <j>  from  eq.  3.4.  The  form  of 
*— which  is  most  useful  in  this  case  is  given  by  eq.  C.5.  As  it 


r  *>0 

is  stated  in  Appendix  C, 


— cosG  -  — 

1  Q>P  =  A _ r 

r  r2  _  0  r  777 

*  —nr  2— cosG  +  1 

Az  A 


(C.5) 


Substituting  this  relationship  into  the  appropriate  equation  of  eqs.  B.6, 


t '  icose 

Sr  £  -  zjcose  +  i 


(B.7) 


Eq.  B.7  can  be  integrated  to  obtain  (j  as  follows: 


C  --  -  — COS0 

^  I  7^-7— - 

J  ~  2— cos 9  +  1 


A 


f  •  k2Jw‘ 


2rAcos0  +  A2) 


Acos0 


f  _ 

/  r2  -  2r Acos0  + 


A2  (B .  8) 
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From  the  C.  R.  C.  Standard  Mathematical  Tables  (ref.  6), 


P  dx  -  1  In 

r  x2 

-  ^  1 

r  * dx 

'  x(cx2  +  bx  +  a)  2a 

cx2  +  bx  +  a 

2a) 

cx2  +  bx  +  a 

Applied  to  the  first  term  on  the  right-hand  side  of  eq.  B.8, 


(B.9) 


^  ^r(r2  -  2rAcos0  +  A2)  2  ^n£r2  -  2rAcos0  +  A2j+  ^cos ®  ^ _  2rAcos0  +  A2 

Substituting  eq.  B.9  into  eq.  B.8,  the  equation  for  becomes 


f  =  i  InL - Z*— - JI+ 

'  r2  -  2rAcos0  +  A2 


f(6) 


(B. 10) 


where  the  additive  f(0)  term  is  necessary  since  the  expression  is  derived 

by  the  integration  of  the  partial  derivative  of  0  with  respect  to  r. 

1 


The  evaluation  of  the  f(0)  term  is  accomplished  by  determining  — 

r 

from  eq.  B.10  and  using  the  resulting  expression  in  the  second  equation 
of  eqs.  B.6.  From  eq.  B.10, 


l  =  -Asin0 _  d  [f(9)j  1 

r  "30  r^  -  2rAcos0  +  A2  d0  r 


Noting  that 

r  =  A(cos0  +  cot^5sin0) 

and  r2  =  A2(cos20  +  2sin0cos0cot^5  +  cot2^sin20), 
r2  -  2rAcos0  +  A2  =  A2sin20csc^3 
so  that 

I#.-!  c5c6s1,AJ  +  I  d  &W]  . 

r  ^ 9  A  r  r  d&  Jr 
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But,  from  eqs.  3.6, 


2 

cscbsirr^jl 


-I 

A 

Consequently, 

d  [fX9)],  = 

d9  U 

1  ~ 

or  f(9)  =  constant  =  ^ 

Therefore,^  may  be  written  as 


t-  2 


In 


r^  +  2rAcos9  +  a2 


1  " 

f  |  In  K 


or,  finally, 


f*'  5  Ink 


L" 


Krz _ 

-  2rAcos9  +  A2 


(B.ll) 


Eq.  B.ll  is  the  scalar  equation  for  the  contours  of  constant  j) 

— * 

(the  curves).  A  more  descriptive  form  of  eq.  B.ll  is  obtained  by 

converting  to  the  Cartesian  coordinate  system.  From  eq.  B.ll, 


^  W _ . 

r2  -  2rAcos0  +  A* 

or  r2  -  2rAcos9  +  A2  =  K^e"2^  r2 
Therefore , 

(1  .  Ke-2 <t>  )r2  .  2rAcos0  +  A2  =  0 


or 


r2- 


2rAcos0 


•v  ~2<P 

(1  -  Ke  r  )  (1  -  Ke  r  ) 


=  0 


(B.12) 
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The  conversion  to  Cartesian  coordinates  is  accomplished  from  eq.  B.12  by 
recognizing  that  r2  *  x2  +  y2  and  x  =  rcosQ. 

Eq.  B.12  becomes 


x2  +  y2 


(1 


2A _ 

-  Xe~2<t 


-x  + 
) 


Completing  the  square  in  x  and  rearranging  terms, 


9  2A  .A2  .9 

A2 

A2 

X  (1  -  +  (1  -  )*■  +  y 

"  (1  -  ffe-20  )2 

(TTCTT 

T  A  T2  .  o  A2 

1  . 

(B.13) 

.  .  jTfce-itf  )  +  y2  «  h  .  £e-20  ) 

(1  -  £e-20  )‘  1 

Eq.  B.13  is  in  the  general  form  of  a  circle  in  the  Cartesian  coordinate 
system.  Thus,  the  circular  nature  of  the  VfS  curves  in  fig.  3-1  is 
accurate.  The  radii  of  the  circles  are  defined  by 


and  the  locations  of  the  centers  are 


(B.14) 


(B.15) 


Note  that  eq.  B.15  indicates  that  the  centers  of  all  of  the  possible 
circles  which  are  described  by  eq.  B.13  are  always  on  the  x  -  axis  of 
the  Cartesian  coordinate  system  (line  OQ  in  fig.  B-l). 
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.  /( cos6  +  cot^  sinOl^sin^Ocsc^ 

R*A"y  (cos26  +  coy»si»2e)2 — **- 


or  <v 

R  =  Asin0cs<^ 


cos0  +  cot^  sin0 
cos20  +  cot^sin20 


(B.18) 


Dividing  eq.  B.18  by  eq.  B.3  yields 


R 


(Td  (cos20  +  cot^  sin20) 


(B.19) 


There  is  an  obvious  similarity  between  eq.  B.19  for  -5—  and  eq.  B.4 
a  CTd 


for 


—— ?.  In  fact,  eq.  B.19  can  be  analyzed  as  eq.  B.4  was  for 

>w 


the 


D 

smallest  values  of  the  ratio  - -  which  can  be  expected  by  assuming  that 

3 

or-  3  x  10  radians  and  setting  0  to  45°.  Eq.  B.19  reduces  to 


<3-3xlo2>  tan^ 


(B.20) 


Fig.  B-3  illustrates  eq.  B.20  for  0  ^ /3  &  5°.  If  - r  -  10  is  considered 

'  <Td 

sufficient  to  satisfy  eq.  B.16,  it  can  be  seen  from  the  figure  that  the 
second  part  of  assumption  "d"  is  acceptable  if  ^3  2  i°43' .  Therefore, 
the  linearity  of  the  curves  of  in  the  neighborhood  of  P  for  any 

value  of  0  corresponding  to  essentially  all  the  possible  values  of  ^ 
is  acceptable;  and  the  second  part  of  assumption  "d"  is  generally 
acceptable . 
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APPENDIX  C 


=  -  -^-csc0sin2/3 
J  r  A  / 


r  2  2 

— CSC  0COS0  -  CSC  t 
A 

CSC  2°<S  -  2jcos0  +  1) 


icOS0  -  i 

7?  -  2^COS0  +  1 


Substituting  eqs.  C.3  into  eq.  C.5  and  combining  terms, 


cos 

Asin^0 


^0  +  slnUcosOcot^  -  1 
csc^S  (cos0  +  cot^5sin0) 


or 


icscOsinV  S—  ~  COt^  COS6 

A  I  cos0  +  cot^sinO 


The  magnitude  of 


given  by  eq.  3.7,  is  restated  here  for 


Using  the  expressions  of  eqs.  C.4  and  C.6  in  eq.  3.7, 


sc@sin 


/ 

— 

/ 1+ 

sin0  - 

COt^  COS0 

COS0  + 

cot^  sin0 

L 

(C.4) 


(C.5) 


(C.6) 

convenience • 


(3.7) 


Expanding  the  expression  under  the  radical  and  combining  terms  yields 
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APPENDIX  D 


This  appendix  presents  the  derivation  of  the  expression  for 


GMOP 


,  eq.  3.15,  from  eq.  3.9.  To  begin,  eq.  3.9  is  restated  for 


convenience . 


GMOP  =  — rj — r  =  sinOcsc/5  (cos0  +  cot^  sinO) 
A  N/3 1  '  ' 


(3.9) 


(MOP 


¥ 


6 


=  -  sinOcot^J  csc^  (cosQ  +  cot^0  sin0)  -  sin^Gcsc^gt 


Expanding  and  manipulating, 


^  (MOP 

¥ 


-  sin0cos0cosi 


sp  csc^5  -  sin^Ocos^  csc^S  -  sin^Ocsc^ 
Factoring  and  using  cos^  =  1  -  sin ^3  , 


a  GMOP 

4* 


=  -  sinOcsc 


£  cosOcos ^  +  sinOcsc^  (2  -  sin^  )J 


or 


^  GMOP 

¥ 


-  sinOcsc 


^  ^cos0cos^  +  2sin0csc JS  -  sinOsin^J 


Noting  that  cos(y5  +  0)  ■  cosOcos^  *  sinOsin^d, 


^  (MOP 

V 


*  -  S 


inOcsc^J^  cos (^3+  0)  +  2sln9c,c/s] 


and,  finally, 


^  GMOP 


¥ 


-  sin0csc3. 


^  sin^cos(^  +  0)  +  2sin0  (3.15) 


